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1. Introduction 

Fixed order calculations in Perturbative QCD (PQCD) can be used to compute inclu- 
sive quantities in strong processes. The precision of these computations is limited by 
our ability to compute complex Feynman graphs. Yet, even if we could perform com- 
putations at arbitrary order, we would not be able to give predictions for exclusive 
quantities. This known fact follows from the presence of collinear and soft diver- 
gences in fixed order calculations with a definite final state. Only by summing over 
different final states these divergences can cancel, thereby allowing the computation 
of certain (i.e. the collinear and infrared insensitive) inclusive quantities. 

The only frameworks in which exclusive quantities can be computed must involve 
the resummation of an infinite class of Feynman graphs. Shower Monte Carlo (SMC) 
programs perform this resummation in the leading logarithmic (LL) approximation. 
Going beyond the LL approximation is a considerably complex task: the initial hard 
process would have to be implemented at the NLO order, and shower development 
would have to be improved to NLL accuracy, in both the collinear and soft structure. 
There arc however simpler directions in which SMC programs could be improved. 
Namely, one could stick to the LL approximation as far as the shower development 
is concerned, but improve the treatment of the hard emission, in such a way that 
inclusive quantities are predicted with a next-to-leading order (NLO) accuracy. Such 
an improvement (I will call it SMC-I-NLO from now on) would clearly merge all the 
best features of the two approaches. One would have exclusive final state generation 
together with the accuracy of NLO calculations. 

In refs. [1, 2, 3] a method and implementation of an SMC+NLO program has been 
given. Such method (referred to as MC@NLO) is based upon a careful elaboration 
of the NLO results, that has to match certain features of the SMC program. The 
approximate SMC implementation of the NLO corrections must be subtracted to the 
exact NLO result in order to avoid double counting. The method of refs. [1, 2, 3] 
does not require modifications of the existing shower Monte Carlo code. This is a 
considerable advantage, because Monte Carlo programs are large and complex, and 
modifying them is a major task. The method has however few drawbacks: 

• The approach is specific to a particular shower Monte Carlo implementation, 
since the form of the NLO result has to be adapted to it. Similarly, one 
must extract the NLO terms already present in the Monte Carlo in order to 
subtract them. Variations in the Monte Carlo implementation (like the new 
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showering variables proposed in refs. [4, 5]) would require computing again 
both the modified NLO expression and the subtraction term. 



• The Monte Carlo may not be fully accurate in the treatment of the soft region. 
In this case the difference between the exact NLO result and its Monte Carlo 
approximation may have left over singularities that need special treatment. 

• The correction coming from the difference between exact NLO and its MC 
approximation may be negative. Thus, negative weighted events may be gen- 
erated. 

In the present work, I propose an SMC+NLO method that improves over the 
method proposed in refs. [1, 2, 3] on the aspects listed above. This method cannot be 
directly implemented using existing SMC programs, since it requires minor modifi- 
cations to the shower development. On the other hand, it is quite simple, and it can 
also be used to perform matrix element corrections in a more consistent way. It is 
not unlikely that it may be generalized to implement a SMC-I-NNLO method, when 
NNLO calculations will become widely available. It would be desirable that new 
SMC implementations would eventually provide facilities to implement the method 
proposed here. 

The strategy followed in the present approach is the following. If the hardest 
(i.e. the largest p^) emission was generated first in the SMC, one would just have 
to correct the first emission at NLO in order to get an SMC-I-NLO implementation. 
Unfortunately, in general this is not the case. A consistent treatment of soft emis- 
sions is essential in SMC programs in order to obtain the correct results for hadron 
multiplicity distributions. Large angle soft radiation from bunches of coUinear par- 
tons interfere destructively, and the resulting coherent emission is described in SMC 
programs as arising from the parent parton of the collinear bunch. The correspond- 
ing angular ordered shower algorithms can therefore generate soft emissions earlier 
in the shower [6], so that very often the hardest emission is not the first. The basic 
result of this work is the construction of a shower algorithm which is fully equiva- 
lent to the coherent angular ordered shower, but has the hardest emission generated 
first. This construction for timelike showers is described in sec. 4. There I show 
that the hardest emission should be generated with a modified Sudakov form factor, 
and subsequent emissions should be generated according to the standard algorithm 
modified by a pt veto, so that they cannot be harder than the (first generated) hard- 
est emission. Large angle, soft coherent radiation from the particles arising from 
the hardest splitting must be added, which will result in the presence of truncated 
showers associated to that pair of particles. 

The NLO accurate hardest emission can be easily obtained, as I show in Sec. 5. 
One simply constructs an NLO Sudakov form factor for the hardest emission in 
analogy with the corresponding expression in the SMC. 
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Space-like initial state showers are discussed in detail in Sec. 6. They are usually 
implemented in the backward evolution formalism [7], i.e. evolving from the hard 
event back to the incoming hadron. The shower order is from hard to soft, as in 
the timelike case. This fact allows a similar treatment of the spacelike and timelike 
cases. Certain specific problems in the spacelike case, having to do with the choice 
of the parton density scheme, are discussed in detail and shown to be harmless. 

The paper is organized as follows. In Sec. 2 I collect the main formulae for 
the angular ordered shower. In Sec. 3 I introduce a formal definition of the Monte 
Carlo shower, that will be used in sec. 4 to construct the equivalent formulation of 
the shower in which the largest emission is generated first. In Sec. 5 the NLO 
implementation of the hardest emission is discussed. In Sec. 6 the spacelike showers 
will be discussed. In Sec. 7, the comparison with other matrix element corrections 
methods will be discussed. In Sec. 8 I give my conclusions and prospects. 

2. The angular ordered shower 
2.1 Kinematics 

Shower Monte Carlo algorithms are an effective way of summing up a large class 
of Feynman graphs in hard QCD collision, including collinear singularities at the 
leading logarithmic level, and also soft singularities in the double logarithmic region. 
The implementation given in the HERWIG Monte Carlo program is through an 
angular ordered shower [6, 8]. We summarize here its general features. We consider 
the splitting process shown in fig. 1. In the limit of small emission angles we have 
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Figure 1: Kinematic variables for a splitting process. 
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Defining t — we then have 



Pt = Viz(l 



z). 



(2.2) 



The probability for such a split is given by 



dP = F{z,t)dzdt 



asipr) 
2n 



Pij{z)—dz 



(2.3) 
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where P is the unregularized Altarelh-Parisi sphtting function. The variable z is 
hmited by the requirement Pt = Vt z(l — z) > and to is the minimum t for 

a branching to take place. In the angular ordered shower, subsequent branching can 
take place only with an angle 9 less than the previous branching angle. In terms of 
t, this implies that the maximum t value for a subsequent splitting in fig. 1 must be 
less than z'^t along line 1, and {1 — z^t along line 2. 

2.2 Coherence 

In the angular ordered shower it is quite possible to have a large pr emission fol- 
lowing smaller pt emissions. In the sequence of two branchings depicted in fig. 2, 
for example, one can have z very near one in the first branching, and z' = 1/2 in 
the second branching, so that there is a region of < for which the second pt 
(according to formula (2.2)) is much larger than the first one. The line (1 — z)E 




Figure 2: Two subsequent branchings. 

would be soft in this case, and (if line 2 is a gluon) the splitting kernel cq. (2.3) is 
singular in this region. Thus it would appear that a highly off-shell line can have soft 
gluon emission singularity. In fact this is not what happens. The soft gluon emission 
in this region is the result of coherent radiation from lines 3 and 4, and these line 
have in fact small virtualities. 

3. The shower 
3.1 Notation 

We now want to set up some notation for discussing the properties of a shower. We 
call Ei{t) the angular ordered shower originating from the initial value t. In order 
to give a more precise meaning to this notion we introduce the basic configuration 
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states^ 



\ki,mi]...]ki,mi) (3.1) 



where ki and rrii are the momenta and quantum numbers of the particles, with the 
normahzation^ 

I 

{ki, mi, k, mi\k[, m[; k'l,, m[,) = 6i^i, 6^{ki - k'^6^.^^>. . (3.2) 

i=l 

A shower is defined as 

oo „ 

S = ^ ^ d^ki...£kiC{ki,mi;...;ki,mi){ki,mi;...;ki,mi\ (3.3) 

1=1 mi. ..mi 

SO that given the infinitesimal cell around a basic state 

= \k[, m[; k'i„ m^,) d% . . . d^k'i, (3.4) 

the product 

S ■ (3.5) 

is the probability to generate a state in the cell d"^. A final state observable 
g{ki,mi; . . . ; ki, mi) will be represented in this notation by 

oo „ 

^ d^ki...d^ki g{ki,mi;...;ki,mi) \ki,mi;k2,m2;...;ki,mi) (3.6) 

1=1 mi. ..mi 

and the shower is a linear functional over the observable, which yields its average 
value S • G. 



3.2 The shower equation 

The angular ordered shower obeys the equation 

S{ti) = A{ti, to) (I| + / ' A{ti, t) F{z, t) S{zH) §( (1 -z)^t) dt dz (3.7) 

"'to 

where (I| stands for the initial parton {k,m\, and A{ti,t2) is the probability that 
no emission takes place with the ordering variable between ti and ^2- The first 

term represents the "no emission" case. When we write the product of two showers 
E>{z'^t) §( (1 — 2;)^ t ) we just mean the obvious thing, i.e. the formal product where the 
product of basic states is interpreted as a tensor product. We observe that eq. (3.7) is 
fully meaningful only if, when solved by iteration, it is possible to reconstruct its final 

^We stress that these states do not represent quantum states. 

^We do not need to symmetrize the states over identical particles. Physical observables will 
always be symmetric. 
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state. Thus, we implicitly assume that the parton energy and direction is specified in 
the shower, and furthermore, that the azimuthal angle in the splitting also appears 
in the integration. We do not write explicitly these dependences in order to have a 
lighter notation. 

Equation (3.7) can be written graphically as 



tj ^ 

~ '^z^S^-^y^ 

% (3.8) 

where the large blobs stands for the shower, thick lines are the Sudakov form factors, 
and the small blob is the splitting function. The dashed lines represent the connec- 
tions of the shower blobs to the vertices, and no factors are attached to them. We 
have omitted to write indices along the lines and vertices of the shower, specifying 
particle species. They are inessential in the following arguments, and therefore we 
carry out our discussion assuming that we re dealing with gluons only. 

3.3 Unitarity 

Shower Monte Carlos guarantee that given the initial ti one and only one configu- 
ration is reached. Thus the observable G with g{ki,mi; . . . ; ki,mi) = 1 for all /, k, 
and m must yield S ■ G = 1. In other words, unitarity must be satisfied: the sum 
of the probabilities of all shower configurations must yield one. Thus, from eq. (3.7) 
we immediately obtain 

l = A{ti,to)+[ A{ti,t)F{z,t)dtdz , (3.9) 

Jtn 



and with the ansatz 



we get the solution 



A(t,to) =e"-^4^(^'*')'='^'^*', (3.11) 



which is the usual expression for the Sudakov form factor. Notice that unitarity works 
only if the expression A{tj,t) F{z,t)dtdz is an exact differential. This is equivalent 
to saying that A{tj,t) is the probability that no emission has taken place between 
tj and t. Under this condition, the usual probabilistic algorithm of shower Monte 
Carlo can be used to generate events, i.e. one picks a random number r between 
and 1, solves the equation A{tj,t) = r for t, and then generates z with a distribution 
proportional to F{z,t). 
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3.4 Vetoed showers 

In order to become familiar with the notation we will now prove that the standard 
vetoing procedure^ is equivalent to the use of a modified Sudakov form factor. We 
will use this well known result in the following sections, and thus include its proof 

for completeness. 

A vetoed shower is obtained by applying a constraint 6{g{z,t)) to the branching 
process. One generates a branching at a scale t, and, if g{z,t) is negative, the 
branching is rejected, and one generates a new branching starting from t. The 
equation for a vetoed shower is thus 



^v{ti) = A{ti, to) (I| + f A{ti, t) F{z, t) e{g{z, t)) Sv{zh) ( (1 -zft) dt dz 

J to 

+ I ' A{ti, t) F{z, t) [1 - e{g{z, t))] Sv{t) dt dz . (3.12) 

Jto 



rti 
'to 

Using the identity 



A(*,.*)^^. (3.13) 

dividing both sides by A{ti,to) and taking a derivative with respect to tj we get 

+ / n^.t,) [1 - e(g(z, t,))] dz . (3.14) 



Defining now 



r{ti, to) = exp ^ ' F{z, t) [1 - e{g{z, t))] dt dz^ , (3.15) 



we have 

dr{ti, t 



dt 

so that eq. (3.14) becomes equivalent to 
d r{ti,to)Sv{ti) _ r{ti,to) 



^ = -r{ti, to) J F{z, ti) [1 - 0{giz, ti))] dz , (3.16) 



^ j F{z, tj) e{g{z, tj)) M^'tj) Sv{ (1 - z)' tj)dz. 



dti A{ti,to) A{ti 

(3.17) 
Defining now 

A'{ti,to) ^^^^ = exp(- r F{z,t)e{g{z,t))dtdz] , (3.18) 



r{ti,to) V Jto 



'See for example ref. [9]. 
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and integrating eq. (3.17) using the initial condition of eq. (3.12) we get 



Sv{tj) = A'{ti, to) (I| + / A'{tj, t) F{z, t) 9{g{z, t)) Sv{zH) §y ( (1 - zf t)dtdz. 

(3.19) 

which is a shower equation with a modified Sudakov form factor which matches 
precisely the splitting vertex. 

In the following, we will encounter showers with Sudakov form factors and split- 
ting vertices containing a theta function that limits the transverse momentum. These 
showers can be simply implemented by vetoing, without the need to modify the form 
factor calculation procedure in a SMC. 

4. The largest emission 

We would like now to transform the angular ordered equation in such a way that the 
hardest emission is generated first. The first step to perform such transformation is 
to prove the following three statements: 

(I) the largest emission in an angular ordered shower always takes place along 
the hardest line in the shower. (The hardest line in the shower is the one that 
starts at the beginning of the shower and always follows the largest z line in 
the splitting vertices.) 

(II) along the hardest line, configurations with non-soft emission before the hardest 
emission are collinear subleading. 

(Ill) along the hardest line, for leading configurations, from the first non-soft emis- 
sion down to the last we have t (i.e. t of order p'^ or smaller). 

The proof of (I) is quite simple. Let us call t, z the emission variables of the largest 
Pt splitting. There is a unique way to walk backward in the shower, up to the initial 
tj. Let us call ti,zi, ...ti,zi all the splitting occurring along this line before t,z, 
according to fig. 3. Consider any splitting tj, Zj before t, z. We must have t < Zjtj. If 



/ / / 

/ / / / 




Figure 3: Hardest occurring at t, z, in a shower initiated at tj. 

Pt > Ptj thentz'^{l — z)^ > tjz'^{l — Zj)'^, which implies {tjz'^)z'^{l — z^ > tjZj{l~Zj)'^, 
and thus z'^{l — z)"^ > (1 — Zj)"^, which immediately implies 1 — zj < 1/4, and thus 
Zj > 3/4 for all j , which proves our lemma. 
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The proof of (II) goes as follows. Suppose there is a non-soft emission at t^, Zf^ 
before the hardest emission. Non-soft means 1 — Zk ~ 1. Also Zk > 1/2, therefore 
from ^/tkZk{l — Zk) < Pt we get < fortiori < t. But by angular 

ordering we must have tk > t. Thus none of the collinear integrals from tk up to t 
has a large logarithmic range, and so these configurations are subleading. 

The proof of (III) goes as follows. The first non-soft emission must be the hardest 
emission itself, or may occur after it, but not earlier because of statement (II). Let 
us call t', z' the values corresponding to the first non-soft emission. We must have 
Vt'z'il — z') < Pt, which implies t' < p^- angular ordering the same must hold 
for all following emissions. 

The equation for the angular ordered shower (3.7) can be formally solved by 
iteration. Using the fact that the largest p^ emission has to be on the hardest line, 
we can write a formal solution where the largest p^ is explicitly present: 




S{tj) = A{tj,to){I\+J2 

l,k=0 

(4.1) 

The z,t vertex has the largest px- 

The graphic expression in fig. 4.1 has the following meaning^: 

• Thick lines are Sudakov form factors; more specifically, they are given by 

A{ti, ti), A{zlti,t2), . . . , A{zl^^,ti+k+i,to) . (4.2) 

• The large solid blobs are S((l — ZiYt"^), 

• The hollow vertex blobs stand for 

2F{zi, U) 9{z, - 1/2) 9{pt - VUz,{l - z,)) . (4.3) 
where we have defined p^ = Viz{l — z). 

• The solid vertex blob (i.e. the largest p^ vertex) is given by the same expression 
as the hollow blob without the theta function. 

• All intermediate Zj, U are integrated. 

The factor of 2 accounts for the fact that the largest p^ can occur on either line of 
the splitting. Larger p^ emissions cannot occur along the 1 — Zi (and 1 — Zi) lines, 
according to our statement (I). 

^In the following we will also denote collectively the whole z,t variable sequence as Zi,ti with 
i = 1, . . . , / + fc + 1, where z,t = 2i+i,ti+i, and = 
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Equation (4.1) cannot be implemented as it stands in a Monte Carlo algorithm: 
the Sudakov form factors do not match with the splitting vertices, because of the 
Pt theta function. We thus perform the following manipulation. We rewrite the 
Sudakov form factors as 

A /-t 2 , \ - ff '^^ dt' f dz'F{t',z') 

A{tiZt,ti+i) = e ^ ^ ' ^ (4.4) 

= e" ^1/2 Hpt-z'{i-z')V¥) ^ ~ ^ dt' dz'2F{t',z')e{z'ii~z')V¥~pT) 

The first factor on the right hand side of eq. (4.4) matches precisely the splitting 
vertex in eq. (4.1). The second factor is a remnant. Because of statements (II) and 
(III) the remnant can always be written as 

^- It:Sdt' Jl^^dz'2Fit',z')e{z'(l-z')VF-pT) _ ^~ Jl^^^dt' J dz'F(t',z')eiz'(l~z')VF-pT) _ 

In fact, assume that k is the first non-soft vertex. Because of statement (II) in all 
previous vertices we can replace Zi 1, and eq. (4.5) holds. According to statement 
(III), starting with the first non-soft vertex the t values are not larger than Pt, so 
that the exponent in eq. (4.5) vanishes because of the theta function, and also in this 
case eq. (4.5) is trivially satisfied. 

We will call Anitijp^) the product of all remnants. From eq. (4.5) we obtain 

An{ti,p^) e-/*o '^*'/^^'^(*''-')^(^'(i-^')v^-*'-) . (4.6) 

Equation (4.6) has an obvious physical interpretation. It is the Sudakov form factor 
for not emitting a particle with transverse momentum larger than px- It can be used 
to generate the largest px emission first, in association with the F{z,t) dzdt factor 
present in eq. (4.1). We can thus rewrite eq. (4.1) in the following form 

S{ti) = A(ti,to){I\+ / dte{z > 1/2) dz AR{ti,p^)2F{z,t) ^{{1 - z)h) 

Jtn 



CO I I I I oo 



1=0 k=0 

Double lines are the modified Sudakov form factors, defined as 

Ay{t.zlt,^,) = e-/i1'^*'A%'^^'2n*'.^')^(PT-.'(i^.')VF) _ ^4_g^ 

The two graphic factors in eq. (4.7) are now showers, that can be easily implemented 
in a Monte Carlo algorithm. In fact they are angular ordered showers with the 
addition of a px veto. The veto needs to be applied only along the hardest line, but 
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since no emissions with transverse momenta larger than are possible in any other 
branches of the showers, we might as well apply the veto everywhere, without the 
need to follow the hardest line. 

Our final result for the shower is 



S{ti) = A{ti, to) (I| + / dt [ dz Aniti, Pt)F{z, t) Sv{{1 - z)\ p^) 

Jto Jo 

xSvT{ti,t,p^)Sv{zH,p^) , (4.9) 

where 

• Sv((l — zytjp-r) and SvizHjPT) are angular ordered px vetoed showers. 

• SvTiti, tyPr) is also an angular ordered p^ vetoed shower, starting at t/, except 
that along the hardest line the minimal scale for showering is t instead of tg. 
We call it the truncated shower. 

Unitarity of eq. (4.9) works as follows. Since both Syr(i/,i,pT) and Sy(i, io) project 
to 1 by unitarity, the unitarity equation is 



'*() 

which is the analogue of eq. (3.9). We have 

i-ti 



A{ti,to)+ I dt I dzAR{ti,p^)F{z,t) , (4.10) 

Jtn Jo 

)f eq. (3.9), 

dt I dzAR{ti,p^)F{z,t) 
) Jo 

= / dp^ r dt f dz5{p^-^tz{l-z))F{z,t)e-^'^i''^^''''^''^''^'^''^^-''^^-^-^ 

J Jto Jo 

J dp'j- 

The last step follows from the fact that the value of the exponential is 1 for p^ ^ oo 
(because the exponent vanishes), and it reduces to A{ti,to) for = 0, since in this 
case the theta function is irrelevant. 

The first emission in eq. (4.9) is given in terms of a Sudakov form factor that 
depends upon z and t, through the p, variable. It can be easily expressed in terms 
of Pt as an independent variable 



to 



^ f\,^ f^'" ^ f dz e{VU - t) , (4.12) 

^ ^0 JVU)/4: Pt Jo 



where now t = p| / z^(l — z)^ is the dependent variable. Observe that the lower 
bound in p, is \^/A. As discussed earlier, we always have an implicit cut-off at this 
value of Pt that prevents the argument of ctg from becoming too small. 
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We can express the first emission Sudakov form factor in terms of a px integral 



Ai?(i/,pT) = exp 



PT 



/Tr/4 r,,/ 

dp'^^F{t\z')e{tj-t') 



(4.13) 



Now A/j(t/,pT) can be used to generate the largest emission using the standard 
SMC technique: one generates a random number < r < 1 and solves the equation 
AR{ti,pT) for Pt- Alternatively, one can define 



Ai^(p^^PT)=exp 



JpT 



2f' 

dp' Fit', z') 

Pt 



(4.14) 



generate the event according to A^, and implement the 6{ti — t') by vetoing. 

The Sudakov form factors in SvTiii-ii-iPT) and Sv'(^;Pt) include the px theta 
function. They depend therefore upon two variables, t and p^. They can be eas- 
ily implemented using the standard Sudakov form factor (i.e. without the theta 
function), and a px veto procedure as in eq. (3.12) with g{z, t) = p^ — \/iz{l — z). 

The truncated shower is a soft shower. It does not appreciably degrade the 
energy entering the hardest emission. This fact is necessary in order to generate the 
hardest emission first. It is also natural: the radiation of the truncated shower is in 
reality coherent radiation from final state particles, and so it should not steal energy 
from the hardest splitting. 

Observe that hard radiation in the truncated shower does not have any logarith- 
mic enhancement (i.e. neither soft nor coUinear). Formally, this works as follows. 
If z is not near 1 we can only have coUinear logarithms. Because of the angular 
ordering and the pr veto we have 



9<9i, E OiZii^ - Zi) < (ziE) 9 z{l - z) < 



z.:Ee 



9<9, < 
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4(1 - Zi) 



(4.15) 



(which forces z > 3/4). In order to have a logarithmic integral we need a large 9i 
range, and this is possible only if z — > 1. Thus no logarithmic enhancement is present 
for hard radiation. 



5. NLO Corrections 

5.1 NLO expansion of the SMC 

The argument given in the previous Section was referring to a single jet being pro- 
duced. The extension to the general case is straightforward. If there are m primary 
partons, one introduces the total Sudakov form factor 

Ai^({t/},PT) = n ^ni^'i^P-) ' (5-1) 

i=l,m 
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where with {ti} wc denote the set of the initial showering variables t], i = 1 . . .m for 
all primary partons. The probability distributions for the hardest emission is given 

by 



da = B{pi . . .pm)d^r 



X 



AR{{ti},0) + AR{{ti},p-r) J2 Fi{z,t)e{t\-t)dzdt 



=l,m 



27r 



(5.2) 



The derivation of equation (5.2) is a straightforward extension of the argument pre- 
sented in Sec. 4. Notice that in each emission term in eq. (5.2), according to eq. (5.1), 
there are Aji{ti,pT-) factors also for the primary partons that are not emitting. These 
factors compensates for the p^ veto that is applied also to their showers. 

The Monte Carlo generation of the event proceeds as follows. One generates 
the hardest px according to Af>{{tj},p^). The i and z values^ are chosen with a 
probability proportional to Fi{z,t) 6{t\ — t), where t — p\/z'^ /{I — z)"^. Once i and 
z are chosen, we construct the two partons originating from i. We let all the m + 1 
partons shower with a px veto, and add a truncated shower, starting from t\ down 
to t, originating from i. 

Equation (5.2) has the following 0{as) expansion 



da = B{pi . . .pm)d^r. 



= B{pi...pm)d^r 



ft} 

^ / Fi{z,t)dzdt+ ^ Fi\ 

=l,m''*o i=l,m 



z,t)e{t\ 



t) dz — dt 
' 27r 



1+ Fi{z,t)+e{t] 



=l,m 



t) dz — dt 



(5.3) 



where the first two terms in the square bracket arise from the expansion of Ar. The 
+ notation means that the singularities in t and z are regulated according to the + 
prescription, i.e. in such a way that 



/ 

J to 



dtdzFi{z,t)+ = . 



(5.4) 



The + prescription guarantees that unitarity is preserved also at the pcrturbative 
level, through the familiar cancellation of real and virtual singularities. In fact, the 
second term in the square bracket of eq. (5.3) is the virtual correction. Observe that 
also in eq. (5.2), by unitarity, the factor in the square bracket integrates to one, so 
that the cancellation mechanism works to all order in the perturbative expansion. 

5.2 Pcirton event generator 

In a "parton" event generator the 0{as) approximation to the radiation process 
would be implemented by associating to an event with radiation, characterized by 
^We assume for simplicity an isotropic distribution in the azimuth (p of the radiated particle. 
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the variables pi . . .pn, z, t, 0, a counter-event with opposite (negative) weight, and n 
body kinematics pi . . .pn- This association (or projection) is the inverse operation 
of what the SMC does when sphtting one particle into two. In the splitting process 
the SMC substitutes the i*^ particle in the ensemble Pi ■ ■ - Pm with two new particles, 
with momenta computed as a function of pi, z, t and 0. After this, some momentum 
reshuffling is needed, since the sum of the four-momenta of the splitting products 
now has a non vanishing invariant mass, and energy-momentum conservation must 
be enforced. 

The SMC may generate radiation from different primary partons in overlapping 
regions. In this case a radiated parton may be associated with more than one of 
the m primary partons. Thus, in general, for a given event with one emission, there 
can be m projections and m counterterms, one for each association of the radiated 
parton with one primary parton. 

NLO calculations performed using the subtraction method often associate to 
the event with radiation few counter-events with no radiation, according to some 
projections that are chosen for convenience^. 

The 0{as) expansion of the SMC given in term of the subtraction method, as 
presented here, was first obtained in ref . [1] , where it plays a fundamental role in the 
computation of the term to subtract from the NLO cross section in order to avoid 
overcounting. In the present work, our aim will be to substitute the SMC NLO 
result with the exact one, and then work our way backward to the (NLO accurate) 
analogue of eq. 5.2. 



5.3 Exact NLO formula 

The exact next-to-leading expression for the cross section can be described schemat- 
ically by the following formula 



da = B{pi . . .pm)d^m + V{pi . . .Pm)d^r 



+ 



R{Pl ■ ■ .Pm+l)d^m+l -^Ci{pi . . .pm+l)d^rn+l'^i 



(5.5) 



which should be read as follows: the Born and virtual terms generate m body events, 
the real emission term generates m + 1 body events, weighted by R, with some 
associated counter-events with weights Cj, and m body kinematics obtained applying 
the projections to the m + l body event. The projection is required to be IR and 
coUinear insensitive, so that if one parton becomes soft the projected configuration is 
obtained from the full configuration by removing the soft parton, and if two partons 
become coUinear the projected configuration is obtained by merging the two coUinear 
partons. Furthermore, the difference between R and ^ Cj is non-singular. 



specific illustration of this method is given in ref. [10], section 4. 
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5.4 SMC+NLO strategies 

In the MC@NLO implementation [1, 2, 3], the Monte Carlo hardest emission is 
corrected to match the NLO calculation. This is done by the following procedure: 

(i) The NLO cross section formula is rewritten using a projection P that coincides 
with the shower Monte Carlo projection. 

(ii) The shower Monte Carlo approximation to the NLO cross section, eq. (5.3) is 
subtracted to the exact NLO cross section (5.5). This subtraction is possible 
since the two formulae, after step (i), have the same projection. 

(iii) The Monte Carlo is run with its standard input, plus a correction obtained by 
adding a contribution where the subtracted NLO correction of step (ii) is used 
as an input to the Monte Carlo shower. 

The advantage of this method is that it does not require to modify the showering code 
in order to be implemented. Observe that the NLO correction in (ii) is always a hard 
correction, since collinear and soft singularities have been removed from the NLO 
formula by subtracting the Monte Carlo NLO approximation. The disadvantages of 
this approach are 

• The NLO calculation has to be tuned to the Monte Carlo, by changing the 
projection P. It becomes therefore Monte Carlo dependent. 

• One has to determine the Monte Carlo NLO approximation, a task that is not 
always simple. 

• The Monte Carlo implementation of the soft limit may be less than perfect. 

Thus, one has to make sure that the remaining IR sensitivity in the difference 
between the exact and the Monte Carlo NLO expression does not have sensible 
consequences. This last problem is in fact, in practice, a very minor problem. 

• The NLO correction needs not be positive. Thus, one is forced to accept 
negative weighted events. 

If the hardest emission in a shower Monte Carlo is generated first, one has the 
possibility of replacing the first emission with the exact one. We now illustrate 
schematically how this can be done. For case of illustration, we assume that the 
Born final state is described by the variables Vi . . .vi, and the final state with a real 
emission is described by Vi . . .Vi,ri,r2,rs, where the variables are associated with 
the radiated parton. We assume that the projection is simply ¥{vi . . .vi, ri, r2, r^} — > 
{vi . . . vi}. The phase space is written as d^m+i — d^v d^r where is the Born 
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phase space, and d^r is Hdri times a suitable Jacobian. We now write the NLO 
exact formula in the following way 



da = B{v)d^y + V{v)d^y + [R{v, r)d^^d(^r - C{v, r)d%d^ 

R{v, r 



[V{v) + {R{v, r) - C{v, r)) d$,P] d^^ + B{v)d^y 



1 + 



B(v) 



(1-P) d<^r 



(5.6) 



Comparing cqs. (5.2) and (5.3), we immediately see that the analogue of eq. (5.2) 
arising from eq. (5.6) is given by 



da = [V{v) + {R{v, r) - C{v, r)) d$^P] d^^ 

+ B{v)d% 



B{v) 



where we have defined 



(5.7) 



(5.8) 



One can implement eq. (5.7) in an SMC+NLO implementation by generating Born 
events with distribution B{vi . . .vi), generating the first emission according to the 
second line of eq. (5.7), and then generating the subsequent emissions as vetoed 
shower. Furthermore, one should associate a truncated vetoed shower from the 
combined emitted parton and the closest (in p^) primary parton. The first term 
in eq. (5.7) can be generated independently, and attached to an ordinary shower, 
since it is formally of higher order in as- With this method, negative weighted 
events could be generated, since this term is not guaranteed to be positive. A better 
procedure would be the following. One defines 



B(v) = B(v)+V(v) 

+ J {R{v,r)-C{v,r))d'^r 



and then implements the hardest emission as 



da = B{v)d^y 



A 



(NLO) 
R 



(0) 



B{v) 



(5.9) 



(5.10) 



Eq. (5.10) overcomes the problem of the negative weights, in the sense that the region 
where B is negative must signal the failure of perturbation theory, since the NLO 
negative terms have overcome the Born term. 

The structure of the countcrtcrm and the projection in NLO calculations is in 
general more involved than in the example illustrated above. However, one can 
separate the real contribution into several term, each one of them singular in a 
particular collinear region''. To each term one can associate a counterterm with a 
^For example, defining = „ i -j-, where Sk is the mass of the pair formed by the fc"^ parton 

with the radiated parton, we have ^ i?, = ii, and each Rk is singular only in the region where the 
emitted parton is collinear to the fc*^ parton, or soft. 
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projection, and choose v and r variables in such a way that the projection leaves the 
V variables unchanged, and the phase space has the factorized form d^^d^r- The 
discussion given above would go through unchanged, except that in the exponent of 
(5.8) and in (5.10) several terms would appear, one for each singular region, instead 
of the single R/B term, in closer analogy with eqs. (5.1) and (5.2). 

5.5 Alternative formulations 

Several variants of eq. (5.10) are possible. First of all, one may replace B with B 
in the square bracket and in the Sudakov form factor, the difference being of NNLO 
order. 

It is not strictly necessary to have the exact expression for R in eqs. (5.8) and 
(5.7). We can choose an R < R, and add back the positive difference R — R with 
standard Monte Carlo methods, provided that the difference is non-singular, i.e. 
that R has the same singularity structure of R. For example, we may want R to 
smoothly vanish outside of the singular regions, in order to avoid the unnecessary 
exponentiation of large NLO corrections unrelated to the singular region. This is 
easily obtained by defining R — R x hipj-) where /i(pt) goes to one as Pt — and 
vanishes more or less rapidly when is away from zero. There will be a left over 
contribution R — R, which is positive and non-singular, and can be added separately. 

5.6 Running coupling 

In standard NLO calculations one usually uses a coupling constant evaluated at a 
fixed renormalization scale, of the order of some characteristic hard scale in the 
process. The SMC's use instead a running scale for the emission vertex, of the order 
of In order for the NLO corrected hardest emission not to spoil the all order 
features of the SMC, it is convenient therefore to perform an analogous choice in 
eq. (5.10). This can be easily done by taking equal to Pt the argument for one 
power of ctg (i.e. the one associated with the next-to-leading order radiation) in the 
expression for R wherever it appears inside the square bracket of eq (5.10) (i.e. also 
in A^^*^^). By doing this one gets the correct LL Sudakov suppression of events 
with a small largest Pt- It is easy to see that this does not spoil the NLO accuracy 
of the result, since it affects a term which is already of next-to-leading order. 

5.7 Emission from fermions 

Following the hardest line in order to find the hardest emission is only necessary 
for gluon emission lines. Since all emissions before the hardest must be soft, if the 
soft particle emitted is not a gluon one has power suppression of the corresponding 
configuration, which can therefore be discarded. Thus, in the case of a fermion 
line, one can simply follow the fermion line. In fact, by not doing so one generates 
configurations that are not even present in the NLO corrections. The example of fig. 4 
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Fermion of different flavour 



Figure 4: If one follows the hardest line instead of the fermion line, a hardest emission 
which is not present at the NLO level may be generated. 

clarifies this issue. If one follows the hardest line instead of the fermion line one can 
have a hardest emission of a fermion of different flavour, which is not present at the 
NLO level. For the same reason, if one reaches a g qq splitting while following a 
gluon line, one should not proceed any further. In fact, if the hardest emission arises 
after this, according to our statement (II) the 1 — z value of the gluon splitting is 
forced to be small, thus yielding a suppressed configuration, since the gluon splitting 
process is not soft-singular. 

6. Initial state showers 
6.1 Kinematics 

We will now study the Monte Carlo equation for an initial state shower. We denote 
with Si{t, x) the backward shower from a parton of type i, momentum fraction x and 
scale variable t. The kinematic is represented in fig. 5. The emissions are ordered in 
increasing angles towards the hard scattering. Also in the case of initial state showers 
the angular ordering implements interference effects in soft gluon emission. In fig. 5, 
for example, in the region where 9^6' and for soft (gluon) emission, particle 4 is 
emitted coherently from particles 1 and 3, and since the charge difference of 1 and 3 
is the charge of 2 it appears as if it was emitted from 2. 

We consider the region where x (and therefore also all z's) is not small, that is 
to say we ignore the small-x problem now. The evolution variable can then be taken 
to be 



t = E^9^ 




so that 




(6.2) 
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Figure 5: Kinematic variables for a spacelike splitting process. 
6.2 Backward evolution 

In the backward evolution formalism [7] the probability of evolving backward from 
t,x to t',x without emitting is given by [8] 

rr/ / N fix,t') ^ , f(x,t')A(t) 

where / is the parton density. It is thus the same as in the timelike case, except 
for the factor f{x,t')/f{x,t), which accounts for the different probability to find a 
parton with a fraction x of the incoming momentum at the scale t' instead of the 
scale t. As t' — > to the probability Il{t,t',x) goes to its minimum, and for f — > i it 
goes to 1. In order for the formalism to work 11 must also be a monotonic function of 
t'. It is interesting to see how this works, by taking the derivative of Il{t, t' , x) with 
respect to t' 



, d , , Iiit,t',x) f , 



t'F(t', z)f(x, (') + ^P(xlz)f(z, t')-e(z - x) 



(6.4) 



where we have used the Altarelli-Parisi equation for /. Since 

we see that the first term on the right hand side of eq. 6.4 is exactly what is needed to 
turn the regularized splitting kernel P{x/ z) in the second term into the unregularized 
one (thus leaving a positive definite expression in the square bracket) provided the 
argument of ctg in the evolution equation matches the one used in F{t,z). Notice 
that we have also assumed that the evolution variable of the Altarelli Parisi equation 
is the same variable used in the shower evolution. With these assumptions we thus 
have 

^n(t, t\ x) = Ii{t, t', x)—— / -F{t', x/z)f{z, t') . (6.5) 
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From eq. (6.5) we derive the alternative form of the no-branching probabihty in 
backward evolution 



U{t, t', x) = exp I - J^' dt" J ^ F{t", 



x/z) 



f{z,t") 



(6.6) 



which is the form proposed in the paper where backward evolution was first intro- 
duced [7]. 

The spacelike shower equation for backward evolution is given by 



k{t,x) = Ui{t,to,x) + 



(6.7) 



i dt' dz dy Ui{t, t', x) Fij{t', z)S{x - zy)Sj{t', y) fp-(§;(,,-)(i'(l - z)\ y{l - z)) . 

to Ji{X)'(' ) 

The notation S stands for the spacelike shower, and stands for the other parton 
in the splitting ij (i.e., for example l{qq) = g, l{gq) = q, etc.). Observe that the 
starting value of the ordering variable after the splitting is exactly t', which is consis- 
tent with our choice of the ordering variable for backward evolution, eq. (6.1). The 
initial condition for the final state shower is instead consistent with our conventions 
for timelike showers. 

It is instructive to see how unitarity arises in eq. (6.7). One has 

Mx, t) = Ai{t, to) Mx, to) + [ dt' dz dy Ai{t, t') Fij{t', z)5{x - zy) fj{y, t') . (6.8) 

Jto 

Eq. (6.8), together with the definition of the Sudakov form factors* 

i-t J J./ fi 



Aq{t,to) = exp 
Ag{t,to) = exp 



(6.9) 



yields the regularized AltareUi-Parisi equations 



d_ 
di 



Mx,t) 



^ dz 



Fij{t,z)fj{x/z,t) 



(6.10) 



6.3 The SMC parton densities 



Eq. (6.10) is a formulation of the Altarelli Parisi equation which is appropriate for 
angular ordered showers [8, 11]. It differs from the traditional formulation in the 

^The factor of 1/2 in front of Pgg in eq. (6.9) is appropriate for the exclusive splitting functions 

that appear in the Sudakov form factors. It is the symmetry factor associated with the identical 
particles in the final state of the gluon splitting process. For the same reason only the Pgg function 
appears in Aq. In fact, one might replace Pgg [Pgg + Pgg]/2 in A^, and Pgg — > [Pgg + Pqgj/^ in 
Ag, to have more symmetric expressions. 
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definition of the momentum variables, and also in the z dependent argument of a^- 
The corresponding parton density are equivalent to the MS parton densities 

yMS coUinear leading order. They however differ by dominant double log terms 
in the threshold region. In order to compare them, we write the expression of the 
Drell-Yan cross section that they yield. In the HERWIG scheme, the initial evolution 
angle for Drell-Yan is of order 1, and the energy of order Q (the square root of the 
lepton pair invariant mass), so that the initial t value is Q. Using the HW parton 
densities and the Born term for the partonic cross section, we get schematically, for 
the A?"*^ moment of the cross section, 

^^^(Q) = ^^M/r(Q)]' , (6.11) 

where /at stands for the A?"*^ moment of the parton density. On the other hand, in 
the MS scheme we have 

^r(^3) = ^^Vr'(/^)exp|-^ I'dz'^ f %aM') \ , (6.12) 

which also includes double log resummation in the large limit [11]. Equations 
(6.12) and (6.11) have the same dependence 

Q"4r2<(Q) = <iQ)^ f dz'^!^a,{Q\l - zf) , (6.13) 
oQ^ -K Jo 1-z 

which can be derived directly from eq. (6.12) for the MS scheme, and from eqs. (6.11) 
and (6.10) for the HERWIG scheme. This suggests the following identification 

fnt) = /r(/^)exp (-^ f'dz'-^ ^asiA . (6.14) 

n Jo I- z Jt{i-zy r J 

Thus, it is possible to define the HW parton densities in such a way that they resum 
all Sudakov double logs arising in initial state radiation. It is not possible, however, 
to define them to include all double logs. For example, in the DIS process there are 
double logs arising from the final state jet, due to the fact that as x ^ 1 the final 
state jet mass is forced to become small, thus inducing further Sudakov suppression. 
This leads to the well known fact that it is not possible to absorbs soft double logs 
universally in the parton densities. 

If the final state jet is described by an angular ordered shower, the soft gluon 
structure of the final state is also correctly described in the DIS case. However, in 
order to get an accurate total rate in the threshold region, one has to correct the 
parton cross section with the Sudakov term associated to the final state jet [11]. 

When implementing an SMC+NLO program with hadrons in the initial state, 
one should remember that, in order to have NLO accuracy, NLO structure functions 
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should be used. On the other hand, we have seen that for consistency a HERWIG 
scheme should be used for the structure functions. When discussing NLO corrections 
we will have to pay attention to this problem, and we will see in the following that 
it may be overcome in several ways. 

6.4 Hardest emission 

In the spacelike, backward evolution shower it is obvious that the hardest emission 
is to be found along the spacelike line. An arbitrary number of soft emissions can 
occur before the hardest one. They can be implemented (as in the timelike case) 
as truncated, vetoed backward showers. After the hardest emission the backward 
shower will continue with vetoed emissions. 

Having singled out the hardest emission, one has an alternative formulation of 
the angular ordered shower, in which the hardest emission is generated first. It would 
be generated with the HERWIG parton densities, and a transverse momentum Su- 
dakov form factor. Thus, it would include some NLO corrections, more precisely 
the threshold corrections arising from initial state radiation. As discussed above, it 
would not include other threshold corrections, like those arising in DIS from final 
state jets. One can now correct the hardest emission, using the exact matrix ele- 
ments and NLO structure functions. This will guarantee NLO accuracy for IR and 
collinear safe quantities. On the other hand, the rest of the initial state shower can 
be implemented using the LL HERWIG structure functions. This does not affect IR 
safe quantities, and on the other hand guarantees consistency in the treatment of 
soft radiation from the initial state line. 

The equation for the angular ordered initial state shower can be formally solved 
by iteration. The hardest emission can be singled out as in eq. (4.1), yielding 



oo 



t, \Zi,t, \Zi,trZ,t \Zi,t, \ Zk,ti t^ 



§(t,,x) = H(t,,to,a;)(I| + V ' . . . _O_0_q_ q. 

ij^o y' y y^ y^ 

(6.15) 

We will also use the notation Zi/ti/yi for i = 1, Z + /c + 1 to denote collectively the 
Zi/ti/iji, z/t/y and Zi/ii/yi sequences. The graphic expression in eq. (6.15) has the 
following meaning. 

• Thick lines are Sudakov form factors; more specifically, they are given by 

H(t/, ti, x), H(ti, t2, l/i), • • • , H(tfc+i+i, to, Uk+i+i) ■ (6.16) 

• Each y value equals the y value to its left divided by z. Thus yi = x/zi, 
1/2 = yi/z2, etc. 
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• The solid blobs are §(y|(l - ZifU 

• The hollow blobs stand for 



j{Ziyi,ti) 



(6.17) 



• The solid vertex blob (i.e. the largest pt vertex) is given by the same expression 
as the hollow blob without the theta function. 

• All intermediate Zi, U are integrated. 

As in the timelike case, cq (6.15) docs not have a simple Monte Carlo implementation, 
since no-branching factors do not match the branching factors. We can however 
introduce modified form factors, so that they match the splitting probability, and 
collect the leftover. More precisely, using the expression of eq. (6.6) for 11, we rewrite 
it as 



U(t,t',x) ^Uv(t,t',x)xexpl^- J^^ dt" j ^F{t",. 



fix/z,t" 
fix 



^^(V^x(l-^)-p,)} , 
(6.18) 



with 



Uy{t,t',x)=exp!^-J^dt" J ^F{t",z)l^0^e{p,-VFx{l-z))'^ . (6.19) 

Now Hy matches exactly the emission vertex. It can be easily implemented by 
generating the next evolution scale with the usual form factor, supplemented by a 
Pt veto procedure. The remaining factors should be collected together, yielding 



exp 



l+k+l 

-E 

i=0 



i+1 



dt ^ Fit, -')-P^ 



(6.20) 



The Ilrem term is easily shown to be equivalent to 

Jix/z,t) 



Unit I, to, x) = exp 



f'dt- Fit, z)Ml2^ei^txil -z)- Pt) 

J to ^ j[X,t) 



(6.21) 

In fact, the theta function in eq. (6.20) prevents the appearance of any large (soft or 
coUinear) log when t ^ p^. Now, if there is an i for which Zi is not near one, for all 
j > i we have ti < Pt- In fact, at the i*^ vertex we must have 



tiViZiil - Zi) <pT^ti<PT . 



(6.22) 



But for all tj with j > i we have tj < ti because of angular ordering, and so also 
tj ^ Pt- It follows that in eq. (6.20), for all terms in the sum that are not suppressed 
by the 9 function all previous Zi values are near 1, and so the corresponding value 
is near x. All terms in the sum have then the same integrands, and the integration 
ranges can be joined together to yield eq. (6.21). 
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6.5 NLO Correction 



The SMC+NLO correction to the hardest emission is obtained as in the timehke 
case, and in fact it has the same expression given in eqs. (5.7) and (5.8), provided 
one includes in the definition of B, V, C and R the appropriate structure function 
factor. By doing so, in the colhnear regions the R/B expression entering in the 
Sudakov form factor (5.7) matches the expression in cq. (6.21). 

As in the case of the timelike shower, one can use as the argument of one 
power of ctg in the expression for R appearing in the square bracket of eq. (5.10). In 
the same instances one can also replace the NLO parton densities with the HERWIG 
ones, at the price of a NNLO unknown effect, and with the advantage of a treatment 
of the all order double log region which is consistent with the original SMC. Thus, it 
is possible to perform the matching in a fully consistent way. Only practical studies 
of real implementations will tell whether all these subtleties are of any importance 
at all. 

7. Comparison with other methods 

I will now discuss comparisons and relations with few popular methods used to 
perform some sort of NLO improvement of the SMC. 

7.1 Matrix element corrections 

With this method one corrects the large angle emissions. Thus, large angle gluon 
emission become accurate at the NLO level, but the full integral of emissions at any 
angle has still only LO accuracy. A discussion of this method appropriate to SMC 
with interfering gluons has been introduced in ref. [9]. One generates the shower as 
usual, and applies a correction factor corresponding to the exact tree level expression 
divided by the SMC approximate one only if the emission is the hardest that has 
taken place so far. In order for this method to work, one has to make sure that the 
SMC fills all the phase space for gluon emission, in order not to miss any region. 

The method proposed in the present work can achieve the same purpose, just by 
using B instead of B in eq. (5.10) so that one just needs to know the matrix elements 
for real emission. It has the advantage that one needs not to worry about filling the 
phase space for gluon emission. 

7.2 Tree level initial state with radiation 

In the SMC generation of e^e~ hadrons events one usually starts with a qq event 
as the hard event. When three jet events are studied, one can start with a qqg 
event, requiring a cut on some final state variable (like thrust) in order to stay away 
from the two-jet region. In SMC including gluon interference, one has then to set 
appropriate initial showering angle for each final state line. This poses no problems 
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when the initial state hnes form large angles. It may cause problems if some partons 
form relatively small angles. To illustrate the problem, we consider the e+e~ — qqg 
example. We draw the colour lines for this process in the planar approximation, as 
shown in fig. 6. 




Figure 6: Radiation from the qqg system in the planar approximation. 

The gluon has a colour and anticolour line associated with it. Thus its emission 
angle can be controlled either by 9i or by 62- On the other hand, in the angular 
ordered formalism one generates a single shower starting from a given parton, with 
a given initial angle. The HERWIG program chooses the initial condition for the 
evolution of the gluon jet by picking 61 or 62 with equal probability. In both cases 
the radiation probability of the gluon line is given by the gluon coupling Ca- This 
works properly for the next soft emission in the shower. In fact, let us examine 
the soft emission pattern of the qqg system in the large limit. It is given by 
the incoherent sum of the emissions from the two antennas 01,02 and &i,62- The 
h antenna only radiates inside the small angle (i.e. ^1) cone in the qg direction. 
The a antenna fills the whole solid angle. The colour charge of the gluon's colour 
lines 02 and hi is Cp- Near the gluon, for emission angles below 9i both antennas 
give a coUinear enhanced contribution of the same size, and their sum correspond 
to replacing Cp — > 'iCp ~ Ca- In HERWIG, near the gluon (i.e. for angles in the 
9i cone) one has in fact radiation with charge Ca- Near the gluon, but for emission 
angles larger than 6'i, only 50% of the events radiate with Ca intensity. The net 
effect is radiation of the qg system with a. Cp intensity in the large angle region. 
This is what we expect: the a antenna radiates with intensity Cp-, and the h antenna 
does not radiate for angles larger than 9i. This clever scheme, however, does not 
work for further emissions. In fact, lack of soft radiation around the qg system 
should be suppressed by a Sudakov form factor. In the HERWIG scheme we instead 
have just a 50% suppression. In fact, in 50% of the events the qg system does not 
radiate for angles greater than 9i. This means that for pseudorapidity below the 9i 
angular region the qg system would be scarcely populated by soft radiation. It is 
not clear whether this fact may cause observable problems. At the end of the shower 
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the coloured parton associated with the gluon colour line will have a colour partner 
originating from soft radiation from the a antenna. According to the argument given 
above this partner may not be close enough in rapidity to form a low mass colour 
cluster. 

The method proposed in this paper suggests instead to assign an initial showering 
angle of 9i to both the quark and the gluon. One then adds a truncated shower from 
82 to 9i associated with the qg system. It is clear that this procedure gives the right 
amount of soft radiation colour connected to the outgoing gluon line. 

7.3 The CKKW formalism 

The CKKW formalism [12] is a method to merge perturbative angular ordered show- 
ers with exact matrix elements. The method requires a cutoff on the cluster mass, 
which is used to separate the application of exact and approximate matrix elements. 

There are several analogies among concepts that have arisen in the present work 
and in the CKKW paper. In particular, the need to introduce vetoed showers arises 
also there. 

Here I would like to focus on one concept that has no parallel in the CKKW work: 
the concept of a truncated shower. In fact, CKKW avoid the need of a truncated 
shower using a trick^ that I will illustrate here referring to the e^e~ example of fig. 6. 

In the CKKW method one limits the gluon radiation to 6*1, but lets the quark 
radiate at all angles. In the truncated shower scheme one requires that both the 
gluon and the quark radiate starting from the angle 61, and then adds a truncated 
shower of the qg system for larger angles. Since the colour charge of the qg system 
is equal to the charge of the quark, the radiation of the truncated shower combines 
with the radiation of the quark to yield quark radiation in the whole solid angle, 
which is the CKKW method. This also works for gluon splitting: one lets the harder 
gluon radiate at all angles. The g ^ qq has to be treated differently. Here one should 
let both quark and antiquark radiate at large angle, so that one reconstructs the Ca 
charge of the truncated shower by summing up two Cp^s. 

As far as particle flow and multiphcities, the CKKW method is equivalent to 
the truncated shower method. It is however different in the colour connection that 
it assigns to final state partons. The colour connection of a soft emission in the 
truncated shower is the same as if the emission came from the parent parton. This 
is what it should be, since it is coherent radiation. In the CKKW method colour 
connections are not respected, and this may cause problems with large mass colour 
connected parton clusters. 

8. Conclusions 

In the present paper I presented a method for the inclusion of NLO corrections in 
^See Sec. 3.2 of [12]. 
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Monte Carlo programs with angular ordering, that avoid the problem of negative 
weighted event. The method may be applied also for the less ambitious goal of 
correcting the large angle emissions in a SMC. Also in this contexts it has some 
advantages. The implementation of this method requires minor modifications to the 
showering scheme in SMC programs. One modification is the inclusion of pt vetoing 
in the shower. This is common practice in SMC, and should not present particular 
problems. Truncated showers instead have never appeared in SMC implementations. 
They play a central role in the present work, since it does not seem possible to imple- 
ment the soft radiation of a collinear bunch of partons without truncated showers. It 
is therefore desirable that in the context of rewriting or developing new Monte Carlo 
programs will provide facilities to implement truncated showers. Besides opening 
the possibility of SMC+NLO implementations for several processes, this would also 
allow to improve matrix element corrections schemes with a limited effort. 

The method suggested here also hints at the possibility of SMC+NNLO imple- 
mentations. One should isolate the two hardest emission in a SMC event, consider 
again px vetoed showers at the smaller of the two hardest emissions, and add trun- 
cated shower associated with two or three partons arising from the hardest splittings. 
Although it may be premature to discuss this problem in details now, it may be con- 
venient to consider the possibility of implementing truncated shower from more than 
two lines when building SMC programs. 

The method discussed in the present paper has been developed in the context 
of angular ordered showers, in particular as implemented in the HERWIG SMC. 
The aim of this work was to provide an approach along the hnes of ref. [1, 2, 3], 
improved over some substantial aspects. Different approaches to the problem of 
NLO improvement of SMC programs have been presented in refs. [13, 14, 15, 16]. 
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